Chapter 6 in Progress in Ferromagnetic Research 

Editor: V.N. Murray (Nova Science Publishers, New York, 2005) 



Ferromagnetic phase transition in spinor Bose gases 

Qiang Gu 

Institut fur Laser-Physik, Universitat Hamburg, 
Lumper Chaussee 149, 22761 Hamburg, Germany 



Abstract 

The achievement in cooling alkali atomic gases, such as 87 Rb, 23 Na and 7 Li, to quan- 
tum degeneracy opens up a way to study magnetism in spinor bosons, because these 
constituent atoms usually have a hyperfine spin degree of freedom. This article reviews 
several basic problems related to the ferromagnetic phase transition in spinor atomic 
Bose gases from a theoretical perspective. After a brief discussion on various possible 
origins of the ferromagnetic interaction, the phase diagram of the ferromagnetically 
coupled spinor bosons is investigated. It is found that the ferromagnetic transition 
occurs always above Bose-Einstein condensation and the Bose condensate is fully po- 
larized. The low-lying collective excitations of the spinor condensate, including spin 
and density modes, are discussed. The spectrum of the density mode is of the Bogli- 
ubov form and the spin wave spectrum has a /c 2 -formed dispersion relation at long 
wavelengths. The spin-wave stiffness coefficient contains contributions from both the 
"normal" and the condensed part of the gas. 



1 Introduction 

Itinerant ferromagnetism of Fermi (electron) gases has long been a research topic in solid 
state physics [1]. This phenomenon has already been well understood today, although not 
everything is clear. Models based on a picture of almost free electrons is fairly useful in 
describing the ferromagnetism. It is well-known that a free electron gas exhibits paramag- 
netism (called Pauli paramagnetism), owing to the existence of the Fermi surface. If there 
is a magnetization density M in the gas, the Fermi surfaces for spin-up and spin-down 
electrons are split and consequently the band energy would be increased. When an effec- 
tive ferromagnetic (FM) exchange I s is present, electron gases can exhibit ferromagnetism. 
Within the framework of the Hartree-Fock-Stoner theory, I s results in a negative molecular 
field energy when M is finite. As long as I s is large enough, the value of the molecular 
field energy becomes larger than that of the increase of the band energy caused by M. In 



1 



this case, a FM ground state is energetically favored. The critical value of I s is called the 
Stoner threshold. 

Relatively, the magnetism of Bose gases was less studied in history. The reason why it 
was so may be partially owing to the fact that liquid 4 He, the most prototypical Bose system 
being studied in earlier years, is a system of scalar particles and does not display magnetism 
at all. The first example of Bose gases with internal degrees of freedom should be the cold 
atomic hydrogen which received much attention in 1980s. The hydrogen has hyperfme 
spin F and thus it had ever attracted some theoretical interest in studying its magnetic 
properties. However the original purpose in surveying this atomic gas is to explore gaseous 
Bose-Einstein condensation (BEC), which is the intrinsic phase transition in the Bose gas. 
Since the ultimate goal of BEC was not yet attained at that time, efforts towards other 
aspects, e.g. the magnetism, were quite limited. 

Great changes took place in 1995, when BEC was experimentally realized in alkali 
atomic gases, such as 87 Rb [2], 23 Na [3] and later in 7 Li [4]. Since then, research works 
on the physics of ultracold atomic gases have grown explosively in the communities of 
atomic physics, quantum optics and many-body physics [5,6]. Meantime, research interest 
has already gone beyond BEC itself. Alkali atoms also have hyperfme spins, as the atomic 
hydrogen does. So magnetism of spinor atomic gases is now among the most active research 
topics in this field. 

Earlier experiments leading to BEC in alkali atoms were performed in magnetic traps. A 
relatively strong external magnetic field H was applied to confine the BEC system. Because 
the atomic spin direction adiabatically followed H, the spin degree of freedom was frozen. 
As a result, the atoms behaved like scalar particles although they carried spin. In 1998, 
Ketterle's group at MIT succeeded in confining the atomic condensate in optical traps [7]. 
The condensate was first produced in an magnetic trap, then transferred into the optical trap 
for further study. In 2001, Barrett et al. realized all-optical formation of BEC in which the 
Bose condensate was created directly in the optical trap [8] . The optical trap is a radiation 
electric field with a intensity maximum in space created by focusing laser beams. If the 
frequency of light is detuned to the red, the energy of a ground-state atom has a spatial 
minimum, and then the atom is confined by the electric field. A far-detuned optical trap 
can confine all the hyperfine states equally. With the external magnetic field sufficiently 
low, the spin degree of freedom of optically trapped atoms remains active. So investigating 
magnetic properties of spinor condensates becomes experimentally possible. 

Theoretical investigation of ferromagnetism in spinor bosons can be traced back to the 
early studies on the atomic hydrogen. Siggia and Ruckenstein proposed that the magnet- 
ically trapped spin-polarized hydrogen could exhibit a coherent ferromagnetism, based on 
a phenomenological model [9]. The spontaneous magnetization is perpendicular to the 
stabilizing magnetic field. Yamada considered an ideal spin-1 Bose gas [10] and argued 
that the spinor bosons could exhibit an intrinsic ferromagnetism associated with BEC. 
He showed that the magnetization M{H) of the ideal spinor Bose gas was finite even 
if H = once BEC took place, which suggested that the system is magnetized sponta- 
neously [10, 11]. He called this phenomenon the Bose-Einstein ferromagnetism. Caramico 
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D'Auria et al. drew the same conclusion for interacting bosons (in which the interaction is 
spin-independent) [12]. These results appropriately reveal that the spinor Bose gas is rather 
apt to be magnetized by an external magnetic field. 

The realization of optical confining of atomic condensate significantly provoked the- 
oretical interest in magnetism of alkali atoms. Soon after the success of the MIT group, 
Ho [13] and Ohmi and Machida [14] studied the spinor nature of the F = 1 atomic con- 
densate. They pointed out that a hyperfine spin-spin exchange interaction coming from 
the s-wave scattering can be the dominant interaction in these gases. The condensate can 
be either ferromagnetic or "polar", depending on the sign, but regardless of the value, of 
the exchange interaction. They also discussed the collective excitations in these two kinds 
of condensates, including density and spin modes. Spectra of the density mode are of the 
Bogliubov form in both cases, while the spin-wave spectra have different dispersion rela- 
tion. 

The reason why the spinor bosons are so sensitive to the external field and the spin- 
dependent interaction is attributed to the following fact: without spin-dependent interac- 
tions, the ground states of spinor bosons are degenerate and both the ferromagnetic and the 
"polar" states are among the ground states [15]. Therefore an infinitesimal external field or 
spin-dependent interaction can lift the degeneracy, leading to a ferromagnetic or polar state. 

In this article, some recent progress on the ferromagnetism in spinor atomic bosons is 
reviewed. In Section 2 we explore the origin of the ferromagnetism. Various mechanisms 
of the feiTomagnetic couplings between atomic bosons are discussed. The phase diagram of 
the ferromagnetically coupled spinor bosons is studied in Section 3. Critical temperatures 
of both BEC and the FM transition are calculated. The Goldstone modes accompanying 
the two phase transitions are discussed in Section 4, with special attention to the spin wave. 
Our discussions are dedicated to the 3-dimensional homogeneous Bose system in the ther- 
modynamic limit. Experimentally, atomic gases are in a different situation. Some of those 
experimental conditions are briefly discussed in the final section. 

Bosons with internal degree of freedom could be formed in Fermi gases through spin- 
triplet Cooper pairing, because a pair of bound fermions as a whole behaves like a boson. 
The triplet Cooper pairing has been observed in superfluid 3 He and in some solid state 
materials. Although the Cooper pair is not identical to a boson because Cooper pairing is 
not local and different pairs are strongly overlapped in space, one can expect that triplet 
Cooper pairs exhibit similar behaviors to those of the spinor bosons in some aspects. In this 
article, ferromagnetism in triplet-Cooper-paired Fermi gases is also concerned. 

2 Spin-dependent interactions in atomic bosons 

Interatomic forces in dilute atomic gases are rather weak. The spin-dependent interactions 
are even much smaller than the spin-independent one. Nevertheless, they are key ingredient 
in understanding magnetism of the spinor bosons. Up to now, the following spin-dependent 
interactions have been theoretically discussed. In some cases, they can be ferromagnetic. 
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2.1 Scattering between different internal states 

The interaction between atoms V a t (r) should be a function of the separation r of the two 
atoms. However, it is impossible to evaluate the r-dependence of the interaction in detail, 
especially at short separations. When r is comparable to the atomic size, the standard Born- 
Oppenheimer separation of the centers of mass and the internal electronic structures of 
involved atoms may break down. Consequently, atoms can not be regarded as point particles 
any longer and V a t(r) is not even definable [6]. Fortunately, r is relatively large in dilute 
atomic gases, typically of order 10 2 nm, at which the Born-Oppenheimer approximation 
should be good enough and V a t (r) could be well defined and approximated by the van der 
Waals interaction. To avoid having to calculate short range interactions in detail, the true 
interatomic interaction V a t(r) may be replaced by an effective interaction [16], 

U(r = r 1 -r 2 ) = ^-^5fa-r 2 ), (1) 

m 

where a is the low-energy s-wave scattering length. The s-wave scattering perfectly de- 
scribes the interactions of dilute alkali atomic gases, with a is of the order of lOOae where 
as is the Bohr radius. 

In case that involved atoms have internal degree of freedom, for example the hyperfme 
spins for alkali atoms, interatomic interactions may give rise to transitions between different 
sub-states. Two atoms in the initial state \a, (3) may be scattered by collisions to the state 
\a',f3'), where a and j3 denote the hyperfme states of the two atoms. In this case the 
scattering is a multi-channel problem and the scattering length could vary according to 
channels. That is to say, the scattering length a a p )a ip> is "spin-dependent", with a large 
number of free parameters. The number of free parameters could be reduced by considering 
particle exchange, time-reversal and especially spin rotational symmetries. Suppose that the 
collision between atoms does not change the hyperfme spin F (this is usually the case for 
sufficiently low energy collisions if the atoms are in the lowest hyperfme multiplet) of the 
individual atoms, the pairwise interaction keeps rotationally invariant in the hyperfme spin 
space. Then the collisions between two atoms with hyperfme spin F\ and F2 only depend 
on the total spin / = F\ + F2 and Ufa — r 2 ) is reduced, in the s-wave limit, to [13] 

Ufa - r 2 ) = — ff(n - r 2 ) Va/P/, (2) 

/ 

where aj is the scattering length for collisions between two atoms with total spin /, and Vf 
is the projection operator which projects the total spin of the pair of atoms into / channel. 
Vf satisfies the condition Yl,f "Pf = !• F° r two bosons with spin F, f takes the values 
/ = 2F, 2F - 2, 2F - 4, 0. The / = 2F - 1, 2F - 3, ... channels are forbidden owing 
to the exchange symmetry of two bosons. Likewise, 

F 1 -F 2 = ^A / P / , (3) 

/ 
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with A/ = [/(/ + 1) - 2F(F + 1)]/2. Combining Eqs. © and ©, one can get the effective 
interaction in term of the hyperfine spin operators. In particular, the effective interaction for 
the spin-1 bosons is given by [13, 14] 

U(n - r 2 ) = [g + g 2 F(r 1 ) • F(r 2 )] «5(n - r 2 ), (4) 

m 

where g = (a + 2a 2 )/3, g 2 = (a 2 - a )/3, and F = (F x ,Fy,F z ) are 3 x 3 spin 
matrices. The 5(2 term is a Heisenberg-like exchange interaction. Usually the difference 
between a 2 and ao is very small, so g 2 is much smaller than go- 
Why could scattering lengths be spin-dependent? To answer this question, we recall 
that atoms have internal electronic structures. The pairwise interaction of two alkali atoms 
depends on the spin configuration of the two valence electrons. For example, the interaction 
has an attractive contribution when two valence electrons are in the singlet state since they 
can occupy the same orbital. Otherwise, this contribution is absent. 

The exchange interaction can be ferro- or antiferro-magnetic, depending on ao is greater 
or less than a 2 . In principle, ao and a 2 can be calculated directly by solving Schrodinger 
equation. But exact solutions are not get-at-able for many electron atoms. It is possible to 
check the sign of g 2 indirectly by studying spin dynamics of the spinor Bose condensate, 
since behaviors of the ferro- and antiferro-magnetic spinor condensates are quite different 
[7, 17, 18]. Both theory and experiment suggested that g 2 < for the gas of F = 1 
87 Rb [7, 19], and g 2 > for F = 1 23 Na [7,20]. The spin exchange is antiferromagnetic 
for the F = 2 85 Rb and 87 Rb atoms [17-19]. 



2.2 Super-exchange process 

Now let us consider a Bose gas moving on a periodic optical lattice. Such a system can be 
approximately described by the boson Hubbard model [21,22], 

H H ubbard = - y^tq a \ a i + 2 U °^2 ni ( ni ~ 1 ) - ^ 

(ij) i 

Here (ij) labels two nearest neighbor sites on the lattice, Uq is the on-site Hubbard re- 
pulsion. The optical lattice is produced by the interference of laser beams. Two counter- 
propagating laser beams form a standing wave, which acts on the atoms as a periodic po- 
tential, V(r) = ViS\T?kiri, where fcj is the wave vector of the laser. For a given optical 
potential, tuj\ and Uq are readily evaluated [23]. The Hubbard repulsion is proportional to 
the s-wave scattering length. 

As is well-known, the electron Hubbard model is used to describe the Mott insulator in 
condensed matter physics. Provided the Hubbard repulsion is large enough in comparison 
to the hopping matrix tij, the energy band of electrons splits into two branches separated by 
a energy gap. At half filling, the lower Hubbard band is fully filled while the upper one is 
empty and thus electrons are in the Mott insulating state. With the charge degree of freedom 
being frozen, the low-energy behaviors of electrons can be described by an effective spin 
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model. One can derive the effective spin Hamiltonian directly from the Hubbard model 
through perturbation approach, taking the hopping term as the perturbation. To the second 
order of the perturbation, the spin Hamiltonian is described by the Heisenberg model. The 
spin coupling between localized electrons, called the super-exchange, is antiferromagnetic. 

The boson Hubbard model is introduced to account for the superfluid-Mott insulator 
transition in lattice bosons [22, 23]. It is argued that bosons are insulating in the large Uq 
limit, as the correlated electrons do. Following the suggestion of [23], the superfluid-Mott 
insulator transition has already observed experimentally [21]. By analogy, one can expect 
that the spinor boson Hubbard model is reduced to an effective spin Hamiltonian, too. 

Let us first look at a two-component Bose system, which can be mapped into a pseudo- 
spin-^ boson Hubbard model. Yang et al. [24] and Duan et al. [25] have derived the effec- 
tive spin Hamiltonian, 

tf = J]T Si- Sj, (6) 

which is just the Heisenberg model. However, the super-exchange interaction between 
localized bosons is ferromagnetic, J = —4t 2 /Uo, as opposite to electrons. 

For bosons with integer spins, the effective spin Hamiltonian is relatively complicated. 
Imambekov et al. has investigated the spin-1 boson Hubbard mode in detail [26]. To order 
of t 2 /Uq, the spin Hamiltonian consists of two terms, 

H = J^Y, S i •S j + J 2 ^(S i -Sj) 2 , (7) 

(v) <y> 

with Ji = J 2 = —2t 2 /Uo. The first term tends to stabilize a ferromagnetic order, while the 
second term favors a local singlet of the two coupled spins. 

In above derivations we assume that the number of bosons is the same as the number 
of lattice sites, thus each site is occupied by only one particle. But for bosons, any oc- 
cupation is allowed. The system could be a Mott insulator as long as the boson density 
is commensurate. Ref. [26] presented detailed discussions on the spin exchange at vari- 
ous commensurate occupations. Moreover, an antiferromagnetic on-site coupling between 
bosons is considered in Ref. [26]. J2/J1 increases with increasing either the on-site anti- 
ferromagnetic coupling or the occupation number, and the system tends to have a nematic 
rather than a ferromagnetic ground state. 



2.3 Magnetic dipolar interaction 

Corresponding to the (hyperfine) spin degrees of freedom, spinor atomic bosons have a 
magnetic moment m = 7S, which induces the magnetic dipolar interaction between parti- 
cles, 

Umd = -r^M m i • m 2 - 3(mi • r)(m 2 • r)], (8) 
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where /xo is the vacuum permeability. The influence of the magnetic dipolar interaction 
on the properties of spinor Bose-Einstein condensates has recently attracted considerable 
interest [27-30]. Although very weak, it is expected to largely enrich the variety of phe- 
nomena in ultracold Bose gases, especially in the gases of atoms with larger magnetic mo- 
ments, such as europium where m = 7/j,b (^b is the Bohr magneton) and Chromium where 
m = GfiB [31]. 

The magnetic dipolar interaction is anisotropic and its strength depends not only on 
the separation but also on the spin configurations of the two interacting particles. Pu et al. 
studied the magnetism of an ensemble of spinor condensates confined in a lattice [28,29]. 
The "minicondensate" at each site behaves as a localized mesoscropic magnet and interacts 
with each other via the magnetic dipolar interaction. They showed that the ground state is 
ferromagnetic in a chain, and antiferromagnetic in a square lattice. 

However it is rather difficult to treat precisely the dipolar interaction between particles 
in a gas, because of the anisotropic feature and the long-range character of the interac- 
tion. Many efforts have been done to study the ground state of the magnetic (and electric) 
dipolar gases (or fluids). It was reported that a ferromagnetic/ferroelectric nematic order 
is favored under certain conditions [32,33]. Recently, the dipolar ferromagnetism has also 
been predicted in ensembles of randomly distributed nano-particles [34]. 

Let's demonstrate why the magnetic dipolar interaction can result in ferromagnetism 
from a mean-field viewpoint [30]. When all the magnetic moments are arranged parallel, an 
effective magnetic field B = /ioM is created inside the polarized body where M is the mag- 
netization density, owing to the superposition of the intrinsic field of all the magnetic mo- 
ments. The energy density of the effective magnetic field is f m = B 2 /(2fio) = /j,qM 2 /2. 
On the other hand, the magnetic moment of spinor bosons does respond to the inter- 
nal magnetic field, so the spin direction should follow B, which leads to an energy de- 
crease: f c = — M • B = —fioM 2 . Therefore the total free energy density is negative, 
fm + fc = —fJ'oM 2 /2. Contrarily, in the "polar" or "equal spin" states the magnetic mo- 
ments of particles are compensated: M = 0. At the mean-field level, an internal reference 
particle can not sense the moments of other particles and the magnetic free energy is zero. 

It is worth noting that the dipolar ferromagnetism could manifest itself more easily in a 
Bose gas than in a Boltzmann gas or a Fermi gas [30]. The ferromagnetically ordered state 
can be destroyed due to the entropy increase at finite temperatures. So a Boltzmann gas can 
show the dipolar ferromagnetism only at low temperatures comparable to the energy scale 
of the dipolar interaction. Since the Bose condensate has no entropy, the ferromagnetic 
state could survive in the Bose condensed particles below the BEC temperature which is 
relatively high in comparison to the energy scale of the dipolar interaction. As for the 
Fermi gas, it can hardly display the dipolar ferromagnetism, because the dipolar interaction 
is too weak to reach the Stoner threshold for the itinerant ferromagnetism. That's why the 
dipolar interaction plays a less important role than exchange interactions in understanding 
magnetism in solid. The latter is much stronger. 
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3 Phase diagram of ferromagnetic spinor Bose gases 



As suggested in the Weiss molecular field theory (for classical systems) and the Stoner 
theory (for Fermi systems) [1], the ferromagnetic coupling is expected to induce a ferro- 
magnetic transition in Bose gases. On the other side, in a Bose system whose particle 
number is conserved there exists an intrinsic phase transition, Bose-Einstein condensation. 
It is natural to suppose that the FM transition temperature Tp depends on the energy scale 
of the coupling /. Therefore Tp is possibly smaller than T c , the BEC temperature, for small 
values of /, and Tp > T c if I is sufficiently large. Is that true? 

In this section we study how the FM transition and BEC emerge in a spinor Bose system. 
We first calculate the phase transitions by dealing with a simple microscopic model, then 
analyze the predicted phase diagram from a phenomenological point of view. For simplicity, 
we consider a 3D homogeneous spinor Bose gas. 

3.1 A microscopic model 

In Ref. [35] we have investigated the interplay between the FM transition and BEC. The 
calculation starts with a simplified Hamiltonian given by 

H = -* E a *W - Ih E s i • S i- ( 9 ) 

a (ij) 

The particles are treated as being on some kind of "lattice", (ij) denotes two nearest- 
neighboring sites. The first term in the Hamiltonian represents the kinetic energy. In fol- 
lowing calculations it is replaced by the kinetic energy of free particles with mass m*, 
e k = h 2 k 2 /2m*. 

The second term describes the FM coupling, which could be decoupled via the mean- 
field approximation, 

-E 8 *- Si « - E«S*> • Si + S, • (Si) - (Si) ■ (Si)), (10) 

(ij) {ij) 

where Mj = (Sj) serves as the ferromagnetic order parameter, and is chosen to be along 
the z direction, (Si) = (0, 0, Mj) and Mj = (Sf) = (r(a\ Oia) where aj CT annihilates a 
boson with spin quantum number a at site i. We take Mj = M for a homogeneous boson 
gas. Following the Stoner theory for fermion gases, we call H m = ZI^M the molecular 
field and I s = ZIh the Stoner exchange, where Z is the effective "coordination number" 
which is an irrelevant dimensionless parameter of order unity for a gas. Then the effective 
Hamiltonian reads 

H — Nfj, = E( ek ~ A* ~ vHm)nk a , (11) 

kcr 

where \i is the grand canonical chemical potential. 
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The mean-field self-consistent equations consist of two equations, 



n v^ v '™ n V 

ko ko 



n = N/V is the density of particles, V is the volume of the system and N is the total 
number of particles. M is the magnetization density. Let us consider a F = 1 Bose gas 
with the self-consistent equations given by 



1 = n + U 
M = riQ + ti 



h(a) + h(a + b) + h(a + 2b) 

2 2 

Ma) - h{a + 2b) 



(12a) 
(12b) 



where a = — (/i + H m ) / '(fcgT), b = H m /(kBT), the reduced temperature and coupling 
are given by t = ksTm* / ' (2irh 2 n 2 /' i ) and I = Igi^^m* / (2irh 2 ), respectively, n is the 
condensate density, no = uq/u is the condensate fraction, M = M/n is the normalized 
magnetization, and f s (a) is the polylogarithm function defined as 

We note that / s (0) = C( s )> th e Riemann zeta function. 

Based on this simple model, we have shown that an infinitesimal FM coupling can 
induce a FM phase transition at a finite temperature above BEC. As is well-known, the 
BEC critical temperature t c can be determined by calculating the chemical potential [36]: 
a > for t > t c and a — > as t — > t c from above. Assume a FM phase transition 
is induced by / at the transition temperature tp. We first suppose / is very large so that 
tp > t c . Provided that the FM transition is continuous, i.e., 6^0 with t — > tp, the 
mean-field equations become 

1 =3*|,/3 M, (14a) 

2 

1 = 2/4/1 M- (14b) 
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where cif = a(tp). These equations define a relation between tp and /. / is a monotically 
decreasing function of ap. As ap — >• 0, fz(ap) — > C(3/2) ~ 2.612, and fi(ap) « 



y^ir/ap. So for small values of / and ap we have 

ap « 4vrt / 2 , (15) 

where to = l/[3C(3/2)] 2 / 3 is the reduced BEC critical temperature for the Bose gas with 
/ = 0. ap is larger than zero at tp, which means BEC does not yet take place. The Bose 
gas undergoes BEC at a lower temperature. 
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Figure 1: Reduced temperature t vs. reduced FM coupling I phase diagram of spin-1 
Bose gases. The paramagnetic normal phase (PM), ferromagnetic normal phase (FM), and 
coexisting phase of ferromagnetism and Bose condensation (FM+BEC) are indicated. 

By solving the mean-field equations numerically, we can obtain that both the FM transi- 
tion temperature and the BEC critical temperature increase with the FM coupling. Asymp- 
totic analysis at very small I shows that 



5t F 



st c 

— oc 1, 
to to 



(16) 



where Stp = tp — to an d 5t c = t c — to- We should note that the asymptotic analysis is 
not appropriate for this simple model within mean-field approximation, because numerical 
results show that both FM transition and BEC are discontinuous at very small /. Approxi- 
mately, the critical value of / under which the transition becomes discontinuous is 0.35 for 
the FM transition and 0.2 for BEC. For larger couplings, the FM transition is well described 
as being continuous. This point is consistent with the Weiss molecular field theory for clas- 
sical particles, in which the FM transition is continuous [1], because for large couplings, 
the FM transition occurs at a relatively high temperature, when the Bose statistics reduces 
to Boltzmann statistics. 

Figure 1 shows the phase diagram of spin-1 Bose gases. The curves are obtained by 
solving the reduced self-consistent equations supposing both FM transition and BEC are 
continuous. Diamonds and circles denote numerical results of Eqs. (I12at and <!12bb . In 
order to make a comparison with the Fermi gas, Fig. 2 plots the schematic phase diagram 
for both the Bose and Fermi gases together. 



3.2 Phenomenological analysis 

In this subsection, we examine the above-predicted phase diagram on the basis of Ginzburg- 
Landau (GL) phenomenological theory [30]. The advantage of the GL theory is that one 
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Figure 2: Schematic relations between transition temperatures and FM couplings I s . Tp 
and Tp represent the FM transition temperature for Bose and Fermi gases respectively. T c 
and To denote the BEC critical temperature for spinor bosons with and without couplings 
respectively. Iq is the Stoner threshold. 



can describe properties of a system by analyzing the symmetry of the system in a simple 
way, without a detailed knowledge of the microscopic background. It is applicable to the 
case of interacting bosons, in which the phase transitions, both FM transition and BEC, are 
supposed to be of second order. 

To begin with, we derive an appropriate GL free energy density functional that describes 
the coexistence of BEC and ferromagnetism. Generally, such a free energy density consists 
of three different parts [30, 37] : 

/ t (*,*t jM ) = / 6 (*,*t) + /m ( M ) + /c (* > *t >M ) j ( i7) 

corresponding to the Bose condensed phase, the ferromagnetic phase of the normal gas 
and the coupling between the two phases. Here *t = (* T )* = is 
the complex order parameter of the spinor Bose condensate. The condensed Bose gas is 
described by the two-fluid model. We suppose / m (M) only represents the ferromagnetic 
phase of the normal gas, with the order parameter M proportional to the magnetization 
density. 

Following Ginzburg and Pitaevskii [38], the free energy density of an isotropic spin-F 
Bose-Einstein condensate is modelled as 

fb = • V* + a|*| 2 + ^|*| 4 + ^* a %,F ay • F CTV 3y* 7 , (18) 

2m 2 2 

where |\I>| 2 = and repeated sub-indices represent summation taken over all the hy- 

perfine states, a = a'(T — T 6 °) and T c ° is the BEC critical temperature. Both a' and 0o 
are positive parameters, and (3q contains contributions of the spin-independent interactions. 
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The fourth term has S0(3) symmetry, arising from the spin-exchange interactions. j3 s can 
be positive or negative, depending on the exchange interaction antiferromagnetic or ferro- 
magnetic. Due to the [3 S term, the time-reversal symmetry in the Bose condensate should 
be broken spontaneously. 1 Hereinafter we consider the ferromagnetic case (P s < 0). 

As mentioned in the last subsection, a FM transition in normal gas happens above BEC. 
The free energy density for this ferromagnetic phase can be expanded in powers of |M| 2 : 

IMI 2 IMI 4 

f m = c\VM\ 2 + a'(T-T f y-^- + b 1 -^-, (19) 

where a = a'(T—Tf), c, a' and b are positive constants, Tf is the FM transition temperature 
in the normal Bose gas. We suppose that the ferromagnetic normal gas couples linearly to 
the spinor condensate, 

f e = -gM ■ *;F CT7 * 7 , (20) 

with the coupling constant g > 0. This is the simplest coupling term that satisfies the 
physical situation. 

We do not consider fluctuations in this subsection, so all the parameters are supposed to 
be given by their average: {^ a ) = &a and (M) = (0, 0, Mo) where the FM order parameter 
is chosen to be along the z direction for convenience. &t&cr = no is j ust the density of 
condensed bosons. Hence the gradient terms in Eqs. (II 8b and d 19l > can be dropped for the 
homogeneous system. 

Minimizing the total free energy f t with respect to $>' , one gets 

[ a '(T - T c °) - gM a]$ a + /? |*| 2 ^ + 0sK*°i ■ F„ v $ y $ 7 = . (21) 
Then the stable solution reads 

|$ | 2 = |$-i| 2 = 0, |$i| 2 = — ^— (r-T c °-4Mo). (22) 

Po + Ps v a' J 

Only the spin-1 bosons condense, occurring at an enhanced BEC transition temperature, 

T c = T c ° + ^M . (23) 
a' 

Obviously, the magnetization in the normal gas promotes the BEC critical temperature. At 
T = T c , the order parameter of the condensate is zero, and we can derive the value of Mq 
by minimizing / m (Mo) with respect to Mo, 

M = ^[T f -T c ]. (24) 

'Using the notation (*i,* ,*-i) = (<f>ie iei , <j>o, f^-ie^- 1 ), the j3 s term becomes j3 a 4>l[(j>l + cfi-i + 
2<t> 1 ct>-icos(6i+6-i)]+l3 s (<t>i-(l> 2 - 1 ) 2 /2. Since 9i +9-! = vrfor /? s > and 9i + 9-i = for j3 s < 0, 
it therefore determines the ground state of the condensate: it is ferromagnetic for f3 s < and "polar" for 
fa > [13,14]. 
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We now derive the phenomenological relations between the two transition temperatures. 
At very small ferromagnetic couplings, 1^0, both 5T C = T C — T c ° and 5T f = T f - T c ° 
tend to zero. Substituting Eq. (I2"4l into d2"3l . one finds 

5T C = ^T*(5T f -5T c ), 

with T* = (g/a') 2 a'/b. Under the condition that ST f « T* we have 

m) « STf (l - ^) • (25) 

Suppose STf increases linearly with the ferromagnetic coupling /, 5Tf = CI, when / << 
1, the coupling dependence of 5T C is given by the formula 

5T = CI (l - . (26) 

Besides a linear term, ST C also depends on / quadratically. So far, the phase diagram 
predicted in the previous subsection is roughly reproduced. It is noteworthy that the linear 
dependence of 5Tf on / is only an assumption in phenomenological theory. To decide the 
precise relation between 5Tf and I is definitely of theoretical interest and is still an open 
question. 

In analogy to the ferromagnetically coupled spinor bosons, a similar system has been 
found in condensed matters, say, the ferromagnetic superconductors [39,40]. In these mate- 
rials, such as UGe2 [39] and ZrZn2 [40], coexistence of ferromagnetism and superconduc- 
tivity has been observed. The phase diagram, see Fig. 3, indicates that these materials first 
undergo a FM transition, then go into the superconducting state. Considering the strong 
magnetization in these materials, it is suggested that the Cooper pairing should be p-wave 
triplet, thus behaves like spin- 1 bosons in some sense. Phenomenological theory has been 
performed to explain the phase diagram of FM superconductors [41,42]. This subsection is 
just an extension of the phenomenological theory for FM superconductors to spinor bosons. 

4 Spin waves in spinor Bose condensates 

Once a continuous symmetry is spontaneously broken, gapless Goldstone modes are ex- 
pected. Spontaneous magnetization breaks the spin rotational symmetry and the corre- 
sponding Goldstone modes are referred as spin waves. In a ferromagnet, the dispersion re- 
lation of the spin wave takes the form lo s = c s k 2 where k is the wave number, and the spin- 
wave stiffness coefficient c s should depend on the ferromagnetic couplings. For example, 
the long wave-length spectrum of spin waves in a Heisenberg ferromagnet is hw ~ | J\k 2 
with the simple cubic structure, where J is the Heisenberg exchange. 

The /c 2 -formed spectrum can be derived qualitatively from the GL free energy density 
functional for FM transitions given by Eq. ( fl9l . Within GL theory, the Goldstone mode is 
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Figure 3: Schematic Phase diagram of the FM superconductor ZrZn 2 . Tp and Tg are the 
FM and superconducting transition temperatures. P is the pressure. 

interpreted as transverse fluctuations on the average value of the order parameter. The full 
order parameter of the FM phase takes the form M = (M) + SM = (5M x ,5M y , M + 
5M Z ) where 5M Z is the longitudinal fluctuation and 5M X , 5M y are transverse components. 
To linear order in 5M, the free energy is is given by f m = + fm< with 

/° = l -aM*+ l -bM^ (27a) 

f m = cV5M z V5M z + ^(a + 3bM$)5M z 5M z 

+cV5M + V5M- + - (a + 6M 2 ) 5M + 5M^ . (27b) 

Here 5M + = 5M X + i5M y and 5M + = 5M X - i5M y . Below FM transition point, 
dfm/dMo = and we have a + 6Afg = 0. Under this condition, the dispersion rela- 
tion of the transverse mode becomes gapless, 

Huj± = ck 2 , (28) 

while the longitudinal mode is gapped, frw z = ck 2 + 6Mq and thus this mode can be neglect 
at low energy. Here c s = c is a phenomenological parameter. 

To proceed, we consider Goldstone modes in a spinor Bose condensate. This prob- 
lem has received much attention recently [13, 14,43-46]. Since both the conservation of 
particle numbers and the spin rotational symmetry are spontaneously broken in the ferro- 
magnetic condensate, an interesting questions arise: how the spin wave manifests itself 
therein. Ho [13] and Ohmi and Machida [14] studied the zero-temperature collective ex- 
citations in spinor condensates based on the Bogliubov approximation using an equation 
of motion approach. As they pointed out, the density, spin and "quadrupolar" spin fluc- 
tuations are related to S^Bi, S^q, and S^-i, respectively, since 5n = y / no(<5*I'i + 5^i), 
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8M- = y/n^d^Q and 8M 2 = 2 y /no8 i &*_ v A Bogliubov spectrum for density fluctuations 
and a k 2 -formed dispersion for spin fluctuations were derived, 



4irh 2 h 2 k 2 

fuv i = \l e l + 2 (go + gs)n e k , tuv = e k = — — , (29) 

V m 2m 

consistent with available theories of BEC and ferromagnetism. Ueda obtained the same 
results through a many -body mean-field approach [43]. These results hold for a dilute 
weak-interacting atomic gas, in which the condensate fraction is almost equal to one at 
sufficiently low temperatures. 

In the following, we generalize the equation of motion approach to finite temperature 
cases, especially near the BEC point, making use of GL theory [37]. The full order param- 
eters for a condensate are written as ^ = ($* + 8^\, 8^q, £\E , !_ 1 ). Adopting the results of 
Section 3.2 that only F = 1 bosons condense, the free energy density for a ferromagnetic 
condensate reads 

fb = fb+fb, (30a) 
fS = a^ + ^o + A)*! 4 , (30b) 

fb = 7^V,5*;V^ CT + [a + (/3o + /3 s )$i 2 ]^;^ (T 

+ ft)*i 2 (**i + 5*i) 2 - 2/3 s $ 1 2 «5#V*_i • (30c) 

Here (5 = (3q + (3 S . The equations of motion are obtained from f£ according to ihdtS^ a = 
{dfl)/{d5^%). We have 

m {snj = ^ + a + ^ UU + ^ Ui + m) (31) 



\<5tf _i J \(e k + a + /3$i 2 - 2/? s $i 2 )<5* _i ' 
It is easy to get the frequency of 8^\: 



fiwi = V ( e * + a + ^i 2 ) 2 + 2 ^i 2 ( £ fc + « + $ i 2 )- (33) 
Under BEC, /° satisfies the relation df®/d$i = , which yields 

a + /^i 2 = 0. (34) 

This condition guarantees a gapless Goldstone mode in the Bogliubov form. We can also 
derive frequencies of 8^o and <5\l/_i: 

huj = e k + a + f3^ x 2 , Hu-i = e k + a + (j^i 2 - 2/3 s $i 2 . (35) 
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Under condition of Eq. d34t . loq also becomes gapless. Up to now, the results of Refs. [13] 
and [14] are reproduced, except that the meaning of parameters is different. One advantage 
of the phenomenological approach is that it could reveal more evidently the relation between 
the Goldstone mode and the phase transition. 

However, we should note that the above derivation is not self-consistent, because it ne- 
glects the magnetization in the normal part of the gas. According to discussions in Section 
3, the FM transition takes place above BEC. So near the BEC temperature, the magnetiza- 
tion in normal gas is much larger than in the condensate. Apparently, this FM phase in the 
thermal cloud (described by f m ) should be taken into account. In this case, the coupling f c 
between the two phases plays an important role. f c consists of two terms, 

f c = - 5 M $i 2 , (36a) 
f c = - 5 M (^<5*i - S9*_ X S9- X ) 
/2 

--^-g$ 1 (6^5M + + <f* (SM_). (36b) 

The self-consistent solution to phase transitions and Goldstone modes in spinor condensates 
should be acquired by treating the total free energy ft = f m + ft + fc as a whole [37]. The 
obtained results suggested that the Bogliubov mode remains unchanged, but spin waves in 
the thermal cloud and in the condensate are coupled together, 




This equation indicates that the transverse spin fluctuations in both phases are gapped solo. 
But after considering the coupling, the Goldstone theorem is recovered. The spectrum for 
this coupled mode at long wave length is given by 



gM ck* +J|£^ , 2 1 

fl*ii 2 Ml 



hu « « ck z + hje k . (38) 

A/ r 0$1 2 O 11/2 « 



Once again, we obtain a A; 2 -formed spectrum for spin waves. 

Many attempts have been made to evaluate excitation spectra at finite temperatures mi- 
croscopically, on the basis of perturbation theory [44, 45] and generalized self-consistent 
Hartree-Fock theory [46]. The obtained results show that the dispersion relations are in 
the same forms as at zero temperature. The phonon velocity of the Bogliubov mode de- 
creases while the spin-wave stiffness c s increases as the temperature is increasing. But 
these theories are not very much applicable near the BEC temperature, not only because of 
the invalidity of the calculating method itself at the temperature regime under discussion, 
but also because of the neglect of the normal FM phase. 

As a comparison to the ferromagnetic condensate, let us take a glance at the spin waves 
in antiferromagnetic Bose gases in which the Bose-Einstein condensate is in the "polar" 
state: only spin-0 bosons condense. In this case the spin and density waves are described 
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by 8n = ^{S^o + 5^* ), 6M- = y/nHiSV-! + 5V$) and 5M + = 8 Ml [13, 14]. The 



At very small k, the spin-wave spectrum Hu± ~ ^ffhk. The dispersion relation is linear in 
k which coincides with that in a Heisenberg antiferromagnet in which hw ~ | J\k. But, the 
dependence of the spin wave velocities on the couplings are different. 

5 Conclusions and Further discussions 

In summary, some theoretical results concerning ferromagnetic transitions in spinor Bose 
gases are reviewed. The subject covers the origin of the coupling (which induces the tran- 
sition), the transition itself (including the critical point and its relation with BEC), and the 
Goldstone mode (as a result of the transition). A unique feature of the FM transition in 
bosonic systems is that it takes place always above Bose-Einstein condensation, regardless 
of the magnitude of the coupling. The spectrum of the spin wave takes the same form as in a 
conventional ferromagnet, u ~ k 2 , while the spin- wave stiffness coefficient consists of two 
different parts: one coming from the thermal cloud, the other from the Bose condensate, 
which embodies the "two-fluid" feature of the system. 

The above theoretical results are devoted mainly to a homogeneous system in the ther- 
modynamic limit and the actual experimental situations for ultracold atomic gases are not 
well taken into account. In the following we discuss briefly some experimental facts which 
may cause significant differences. 

1) The trapping potential. Atomic gases are experimentally confined by a trapping 
potential which can be approximated as being harmonic. One direct consequence of the 
confinement is that it changes the density of state of the gas and thus the critical behaviors 
[47]. So how the trapping potential affects on the phase diagram of the FM spinor boson 
deserves further study. More recently, Huang et al. [48] studied BEC in trapped F = 1 
spinor bosons with FM couplings, but the FM transition was not considered. 

2) The number of confined atoms. At present, the number of confined atoms is typically 
in the range of 10 6 — 10 8 , far away from the thermodynamic limit. So questions arise: 
do theories derived in the thermodynamic limit hold in this case? Does the concept of 
spontaneous symmetry breaking apply? And to what extent does it apply? The observation 
of BEC suggests that the spontaneous symmetry breaking is still a valid concept, and the 
observed critical points agree with theories quite well [49]. But it is still an open question 
whether it is so for the FM transition. 

3) The spin conservation rule. It is observed that the total spin approximately conserves 
through the evaluation of the Bose-Einstein condensate [7, 17, 18]. From Ref. [17, 18], 
this conservation dominates the spin dynamics and the final states of the system. This is 
in conflict with the picture of spontaneous symmetry breaking which says the total spin 



spectra are 




(39) 
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does not conserve in the ground state. It might be a consequence of the departure from the 
thermodynamic limit. 

Isoshima et al. studied phase transitions in F = 1 spinor bosons under the spin conser- 
vation rule [50]. They argued that the system may have two spatially phase-separated Bose 
condensates with mp = 1 and — 1 respectively. 

Although many work has been done, the research on ferromagnetism in spinor bosons 
is still in a very primary stage. Many questions remain open. Direct comparison between 
theories and experiments is far from being reached. 
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